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Abstract 



The pp annihilation at rest into the (1)7] final state was measured for three different 
target densities: liquid hydrogen, gaseous hydrogen at NTP and at a low pressure 
of 5 mbar. The yield of this reaction in the liquid hydrogen target is smaller than in 
the low-pressure gas target. The branching ratios of the (f)!] channel were calculated 
on the basis of simultaneous analysis of the three data samples. The branching 
ratio for annihilation into 0// from the '^5'i protonium state turns out to be about 
ten times smaller as compared to the one from the ^Pi state. 
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The Okubo-Zweig-Iizuka (OZI) rule § allows one to relate the ratio between the cross sections 
or annihilation frequencies of and oj meson production to the mixing angle O of the vector meson 
nonet: 

Here 6 — O — Oq, where & = 39° from the quadratic Gell-Mann-Okubo mass formula, Oq = 
35.3'' is the ideal mixing angle and / is a kinematic phase space factor. It is known that the 
predictions of the OZI rule are fulfilled in the hadronic interactions at a level of about 10% (for 
a review, see [Q). Recently large OZI rule violation has been found in some pp annihilation 
channels in the experiments with stopped antiprotons at LEAR (CERN). Measurements of the 
Crystal Barrel collaboration show that the ratio of yields for the reactions p + p — > 4>{uj) + 7 
is R{4>^/uj^) = (250 ± 89) • 10~'^ [||. Large apparent violation of the OZI rule was also found 
for the p + p ^ fpi^) + ""^ channel, where R{(t)TT^ / loti^) — (96 ± 15) • 10^^ for annihilation in 
the liquid hydrogen target |^ and R{4>Tr'^ /ojtt^) = (114 ± 24) • 10~^ for annihilation in the gas 
hydrogen target Q . Conservation of P- and C-parities allows the pp (pn reaction only from 
the '^S'l or ^Pi initial states. An interesting dynamical selection rule was found for this channel: 
(j) production is at least 15 times more prominent for annihilation from the spin triplet ^5*1 initial 
state than from the singlet ^ Pi one j|] . 

On the other hand, there are channels of meson production in pp annihilation at rest where 
no OZI rule violation was observed. An example is the reaction 

P + P — > + ?7> (2) 

where the Crystal Barrel collaboration found that the ratio R for annihilation in liquid 

R{4>'q/uiri) = (6.0 ±2.0) • 10"'^ is in agreement with the OZI rule prediction Since reaction (|^) 
proceeds from the same protonium levels as the pp —> reaction, but from isospin 1=0, one 
may expect that the ^77 final state will also be suppressed for annihilation from the ^Pi state. 
Surprisingly enough, the result reported here shows that the situation is opposite: the yield of 
reaction \^ for annihilation from the P-wave is higher than from the S-wave. 
We performed systematic measurements of the channel 

p + p — > K+ + K-+r] (3) 

at rest in a hydrogen target of three densities: liquid, gaseous at NTP and at 5 mbar. Changing 
the target pressure allowed selection of the pp initial state. The fraction of annihilation from the 
S-wave decreases with decreasing target density. According to it varies from 87% to 20% for 
annihilation in liquid and 5 mbar hydrogen, respectively. 

Earlier we reported |Q| measurements of reaction (^) for annihilation in a gas target at NTP, 
though with considerably lower statistics. The measurements of this reaction for the 5 mbar 
hydrogen target are made for the first time. 

The experiment was performed at LEAR (CERN), using the OBELIX spectrometer The 
experimental setup consists of detectors arranged around the Open Axial Field Magnet, which 
provided the magnetic field of 0.5 T along the beam axis. The Time of Flight (TOF) system 
contains two coaxial barrels of plastic scintillators for charged particle identification and the 
trigger. The Jet Drift Chamber (JDC) provides tracking and particle identification by energy loss 
measurement. The High Angular Resolution Gamma Detector (HARGD) is designed to detect 
neutral particles by their decay into 7. Only the information of the TOF and JDC systems was 
used in the present analysis. All data were collected with the trigger requesting two hits in the 
inner barrel of scintillators and two hits in the outer one. 

Two-prong events with the total charge zero and with the length of both tracks more than 
20 cm for each data sample were selected. Kaons were identified by dE/dx measurements. We 
required both particles to be identified as kaons. The number of events in the data samples, the 
number of events that passed quality cuts, and particle identification Nkkx are presented in 
Table 1. 
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Table 1: Summary of the collected statistics. The number of events on the tape, which passed 
quality cuts (two-prong events with the total charge zero and with the length of both tracks more 
than 20 cm) and were identified as kaons are given. 



Target 


On tape (•10'^) 


After quality cuts (•10'') 


Nkkx (-10^) 


Liquid 


10.4 


6.7 


1.81 


Gas at NTP 


6.7 


4.5 


1.06 


Gas at 5 mbar 


9.4 


5.4 


1.71 



The experimental distributions of the X final state for three data samples without 

correction for the acceptance are shown in Fig. 1 (left to right: liquid, gas at NTP and gas at 
5 mbar). The scatter plots of missing mass squared distributions vs Mj^+k- a-re shown in Figs, 
la-c. One can see two horizontal bands from reactions with i:^ and rj. The accumulation of 
events seen on the left side of the plots corresponds to the c^tt" final state. In the right corner 
of the plots there is a spot from pp — > K^K~ reaction and an inclined band from the 7r+7r~X 
final state with wrongly identified pions. To see more clearly the region of these plots around 
the rj band we zoomed them in the centre (Fig.l d,e,f). A blob from the (pr] events is better seen 
in these plots. Figures lg,h,i show the distributions of the missing mass square recoiling against 
the two charged kaons, where the peaks due to reaction (||) are seen. 

To select events from the (prj reaction the invariant mass distribution of two kaons AI^+j^- for 
events with the missing mass around the mass of rj was analysed. It is shown in the middle part 
of Fig. 2, where events with the missing mass interval 0.26 < M^^^^ < 0.34 GeV^/c'^ (centred 
around = 0.3 GeV^/c^) were selected. The left and right parts of Fig. 2 correspond to the 
Mj^+x- distributions for the missing mass intervals below and above the ry mass, where number 
of events with r] is negligible: 0.15 < M^,,^ < 0.23 (left) and 0.37 < M^^^, < 0.45 GeV^/c'^ 
(right). For the sake of brevity we label the missing mass intervals "A", "B" ,"C" in order of 
increasing M^^^,. 

To illustrate the choice of the intervals in Fig. 3 the distributions on M^^^^ obtained by the 
Monte Carlo simulation of reaction (||) and the main background reactions are shown. One 
can see that the interval "B" contains most of the events from the reaction of interest. The 
background reaction 

p + p — > (/) + 7r° + 7r° (4) 

gives the main contribution in interval "A" , whereas the events from the ^Tr^Tr^Tr" final state 
concentrate mainly in the interval "C" . 

The experimental distributions for the interval "B" (Figs. 2b, 2e and 2h) for events with 
missing mass around the mass of rj show clear peaks in the region of the cj) meson. The (j) signal 
for events from neighbouring missing mass intervals "A" and "C" is far less pronounced. 

To describe the peak in M^+k- spectra we used the Breit-Wigner function, corrected for 
phase space, centrifugal factor and convoluted with a gaussian experimental resolution function. 
The background was treated by a smooth third-order polynomial function. 

The (/) peak positions are 1018 ± 1 MeV/c^, 1019 ± 1 MeV/c^ and 1020 ± 1 MeV/c^ for the 
LQ, NTP and 5 mbar data samples, respectively; these values agree well with the PDG value 
for the (j) mass. The widths of the gaussian, which corresponds to the detector resolution, are 
4.3 ± 0.3 MeV/c^, 3.4 ± 0.3 MeV/c^ and 3.7 ± 0.4 MeV/c^ for the LQ, NTP and 5 mbar data 
samples, respectively. They are in agreement with those obtained by the Monte Carlo simulation. 

To evaluate contamination from the background reactions pp (j)X in missing mass interval 
"B" , the number of events with (j) mesons in neighbouring intervals "A" and "C" was determined. 
The background was subtracted assuming smooth behaviour of the squared missing mass distri- 
bution. Monte Carlo simulation of different background channels confirmed this assumption. 

The main background contribution is from the reaction pp — s- ^tt'^tt'^. Using the preliminary 
results for the frequency of this channel obtained by the Crystal Barrel collaboration with the 
liquid target jF(07r°7r°) = (0.88 ±0.15) ■ 10~* one may estimate the number of ^tt^tt" events in 
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Liquid 



Gas at NTP 



Gas at 5mb 




Figure 1: Scatter plots of missing masses vs M^+k- (ajb,c). The same distributions around the 
r] band are shown in (d,e,f); missing mass spectra (g,h,i). The left column corresponds to the 
data at the liquid target, the middle one to the target at NTP and (c,f,i) to the target at 5 mbar. 
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Figure 2: The distributions on Mj^+k- for three intervals of missing mass for three target densi- 
ties: (a,b,c) for the liquid target, (d,e,f) for the gas target at NTP and (g,h,i) for the gas target at 
5 mbar. The left column corresponds to interval "A" (see the text), cohunn (b,e,h) corresponds 
to interval "B", and column (c,f,i) corresponds to interval "C". Dashed lines correspond to the 
(p meson mass. 
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ABC 




Figure 3: The spectra of missing mass recoiling against K^K system for a) pp — *■ 4>fl reaction, 
b) pp ^7r°7r° and c) pp — > (/)7r°7r°7r°, obtained by Monte Carlo. Lines show the interval limits. 



interval "B":7Ve„ =T-Nann-£reg = (0.88±0.15)-10-4.(ii5.9.io6).(0.165±0.010)-10-2 = 17±3. 
It is in good agreement with our estimation of the background, which gives 17 ± 13 events. 

We also verify the contamination from the channels pp (t>X , where X decays into charged 
particles which are not registered by our apparatus. For this purpose, the registration efficiency 
was determined for the reaction pp — + <f>p, when both pions from p decay were missed. It turns 
out to be Ereg = 1-0 • 10"^ With the known frequency oipp (j)p in NTP = (3.4 ± 0.8) • IQ-* 
[|o|, the estimated number of events from this reaction in interval "B" is N^v ~ 0.2 ±0.1. 

The annihilation frequency J- is 



Xann ' ^reg 

where Nann is the number of annihilations in the target and Ereg is the registration efRciency 
that takes into account the correction for all (p and rj decay modes. 

The simulation of reaction (^) shows that the registration efficiency Sreg for annihilation 
from the ^Pi initial state is less than from the '^S'l one. Without information about relative 
abundances of these initial states for annihilation in different targets one could determine only 
an upper and lower limits for the annihilation frequency, corresponding to the initial states ^Pi 
and ^S'l, respectively. 

In Table 2, the number of annihilations Nanm the number of events from the reaction Nev, 
the registration efficiency Ereg and the upper and lower limits for the frequencies of the channel 
pp —> (j)r] for each sample are presented. 

The systematic uncertainties in the evaluation of the annihilation frequencies are about 8%. 
The systematics includes beam monitoring, trigger inefficiencies and the uncertainties due to 
varying the dE/dx kaon selection criteria. 

It is possible to determine the branching ratios B.R. of the reaction pp —^ (prj by using 
additional information from some models of the protonium atoms. Let us divide the number of 
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Table 2: Number of annihilations Nann, number of events N^v, registration efhciency Sreg and 
upper and lower limits of the annihilation frequencies J- for each sample. 



Target 


Nann ' 10" 


Nev 




J^- 10* 


Initial state 


Liquid 


115.9 


68 ±22 


0.97 ±0.02 


0.60 ±0.20 


^^1 








0.58 ±0.02 


1.01 ±0.33 


iPi 


Gas at NTP 


74.5 


146 ± 28 


1.88 ±0.03 


1.04 ±0.20 


^^1 








1.28 ±0.03 


1.53 ±0.29 


'Pi 


Gas at 5 mbar 


85.9 


152 ± 37 


1.68 ±0.03 


1.05 ±0.26 










1.09 ±0.03 


1.62 ±0.40 


'Pi 



(p events N^v in each data sample into two contributions 

N,, = N,,CSi) + Ne.„CPi) (6) 
which are defined, according to Batty 0, as 

N,,('Si) = {l-fp)-^-E{''Sl)-B.R.CSl)-Nann-SregCSl) (7) 
NevCPl) = fp-^-ECPi)-B.R.CPi)-Nann-eregCPl) (8) 

Here fp is the fraction of annihilations from the P-states, factors | and are the statistical 
weights of the ^Si and ^Pi states, respectively. i?('^S'i) and E{^Pi) are the enhancement factors 
which reflect deviations from the pure statistical population of the level. The values of fp and 
enhancement factors are taken from and shown in Table 3. 

The branching ratio B.R. is the probability that the pp system with the definite quantum 
numbers J^'-' of the initial state annihilates into a given final state. Since our measurements 
at three hydrogen target densities give three values of Nev, it is possible to determine the two 
unknown parameters B.R.{^Si (f>ri) and B.R.i^Pi — > (prf) from eqs. 

By using the enhancement factors from for the DRl model and considering for the anni- 
hilation frequency of pp tt^tt^ reaction in liquid the value of (2.8 ± 0.4) ■ 10^'' [0, we obtain 
a different set of values for fp shown in Table 3, which is useful for testing the stability of our 
results. 



Table 3: Fraction of P-wave annihilation fp taken from the analyses of and |11 and the 
enhancement factors E from Wi. 



Target 


_ /p _ 


Enhancement factors 


iy 







E{'Pi) 


Liquid 
Gas at NTP 
Gas at 5 mbar 


0.06 ±0.01 
0.56 ±0.04 
0.84 ±0.03 


0.13 ±0.04 
0.58 ±0.06 
0.80 ±0.06 


0.989 
0.993 
0.985 


0.856 
0.974 
0.999 



The branching ratios obtained from the fit of the events at three densities are: 

B.R.i^Si 077) = (0.76 ± 0.31) • 10""* (9) 
B.R.CPi <pri) = (7-78 ± 1.65) ■ lO""* (10) 

for the parameters choice and 

B.R.i^Si-^ (pri) = (0.82 ± 0.28) • 10"'' (11) 
B.R.CPi 07?) = (7.59 ± 1.55) • 10""* (12) 

for the parameters from [ pT| . Therefore, for both sets of fp the branching ratio ofpp — > (f)?] channel 
for annihilation from the spin singlet ^Pi state is about 10 times higher than the branching ratio 
from the spin triplet '^Si state. An opposite trend was observed for the pp 
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Using the branching ratios obtained above, the annihilation frequencies at different hydrogen 
target densities can be calculated: 



T = {i-u)-1-eCs^)-b.r.Cs,) + u-^-eCp,)-b.r.Cp,). (13) 

The corresponding results are presented in Table 4. 

Table 4: Annihilation frequencies J^-W^ calculated from (^3|) with B.R. from The upper 

numbers correspond to the P-wave fraction calculated with the parameters of 17 1 and the lower 



Reference Liquid 


Gas at NTP 


Gas at 5 mbar 


1 


1 0.71 ±0.07 


1.33 ±0.15 


1.66 ±0.20 




ill 0.67 ±0.07 


1.30 ±0.14 


1.69 ±0.21 



One could see that in spite of the significant difference in the branching ratios from the S- 
and P-wave states the annihilation frequencies increase only by a factor of two from the liquid to 
the 5 mbar sample. It is due to the relatively small statistical weight of the ^Pi state. 

Several experimental studies of the (prj antiproton annihilation channel at rest have been re- 
ported. The results of these measurements are shown in Table 5. The Crystal Barrel collaboration 
measured the pp — > (f>r] channel for annihilation in liquid hydrogen j^] . The absolute value for the 
annihilation frequency was not determined but the ratio to the pp — s- 07r° frequency was given. 
The value in Table 5 was deduced from this ratio and from the annihilation frequency of the 
pp — > (j)TT^ channel measured by the same collaboration [Q. It is in good agreement with our 
result for annihilation in liquid shown in Table 2 and Table 4. 



Table 5: Experimental results of the previous measurements of the pp —^ (j)?] channel at rest. 

— r^rri ■ 



Reference 



T X 10* 



Comments 



I 

ji 



0.66 ±0.19 
0.37 ±0.09 
0.41 ±0.16 
1.04 ±0.33 ±0.05 
0.74 ±0.22 ±0.03 



liquid target 

gas at NTP 
gas at NTP, LX-trigger 
gas at NTP, upper limit 
gas at NTP, lower limit 



The ASTERIX collaboration [|2| measured the (pr] channel for antiproton annihilation in a 
gas target at NTP. Beside it the yield of (H) was measured with the so-called LX-trigger, which 
enriched the data sample with P-wave annihilation events. Our results for the gas target at NTP 
shown in Table 2 and Table 4 disagree with the absolute values of the ASTERIX collaboration. 

The previous measurements of the OBELIX collaboration performed with a gas NTP 
target with another data set and with a lower statistics than in the present paper, are in agreement 
with the results shown in Table 2. 

The main distinctive feature of the pp — > ^77 channel is the increase in the annihilation 
frequency for annihilation in gas at low pressure. In terms of the branching ratios it leads to the 
conclusion that the branching ratio from the spin singlet ^Pi state is 10 times higher than the 
branching ratio from the spin triplet ■^S'l state (see Eqs. (|^ - |l^)). Just the opposite tendency 
was observed for the pp — > ^tt*' channel H, where the branching ratio from the '^Si initial state 
is at least 15 times larger than from the Pi one. 

The unusual spin dependence of the (pn^ production was explained in JTsf under the assump- 
tion that the ss pairs in the nucleon are polarized, as suggested by the experiments on deep 
inelastic lepton scattering (for review of experimental results, see [^]). Following this model, 
the observed strong violation of the OZI rule is only apparent, because it is due to additional 
production of (j) by the connected diagrams that are allowed if the nucleon wave function contains 
a number of ss pairs. The model predicts that the production should be enhanced from the spin 
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triplet initial states, whereas production of the ss pair with total spin S = should be enhanced 
from the spin singlet initial states. Since 77 meson has a substantial ss component, the production 
of the 4>ri final state could be regarded as a production of two ss pairs, one in the spin triplet 
state and the other in the spin singlet state. The extension of the polarized intrinsic strangeness 
model for such case is not straightforward. We note also that the same strong enhancement of r] 
production from the initial spin singlet state was observed in pp —>■ pprj and pn —>■ pnrj reactions 
[15|. An attempt to interpret this effect in the polarized nucleon strangeness model was done in 

In some models opulent production of the (j) meson in pp annihilation at rest is assumed to 
be due to the rescattering diagrams with K*K or pp in the intermediate state (see [ ^) 
and references therein). Calculations of the triangle diagram with K*K intermediate states 
give for the branching ratio of the 4>rj from ^Si initial state B.R.i^Si (prj) = (0.3 ± 0.1) • 10"'', 
a value not far from our experimental results (p[]l^). However, these approaches could predict 
neither the observed spin dependence of the ^tt*^ final state nor the spin dependence of the (firj 
channel. 

An interesting possibility considered in |l^ is that the final state interaction (FSI) of two 
kaons could enhance the (j) production. Indeed, for annihilation at rest the space volume of 
KKrj final state is quite limited, two kaons are created with low relative momenta and could, in 
principle, fuse into (j) due to the FSI. However this model does not explain why the FSI effects 
are stronger for annihilation from the P-wave than from the S-wave. 

An attempt to calculate the (j) yields in pp annihilation at rest in a non-relativistic quark 
model with the ss admixture in the nucleon wave function was done in pc| ]. It was assumed 
that the cf) is produced as a "shake-out" of the nucleon ss component. The calculated branching 
ratios of the t/)?? channel are B.R.{^Si cjir]) = (1.4 - 1.8) • 10"^ and B.R.(^Pi (j)r]) = 
(0.15 — 0.2) • 10^"*. Therefore this "shake-out" mechanism predicts a trend which is exactly 
opposite to that experimentally observed by us. 

The increasing of the (prj yield with the decreasing of the target density arises the question 
of the corresponding behaviour of the pp —^ Lurj channel. This yield was measured only for the 
annihilation in liquid |^ and it is quite high: J-{ujri) — (1.51 ± 0.12) • 10~^. If the arguments of 
the polarised strangeness model are vafid for the 77 production, then the yield of the Y{ujr]) final 
state from the ^Pi state should be higher than from the ^Si final state. It would be interesting 
to verify this prediction experimentally since we have demonstrated that the branching ratio of 
the (j)r] channel significantly increases for the P-wave annihilation. If the ratio between the ^77 
and LOT] channels does not change with the hydrogen target density, the apparent OZI violation 
will be absent for annihilation in low pressure gas as it is for annihilation in liquid. 

In conclusion, the pp annihilation at rest into (prj final state was measured at three different 
target densities: liquid, NTP, 5 mbar. The yield of this reaction for the liquid hydrogen target is 
smaller than for the low pressure gas target. The branching ratios of (prj channel were calculated on 
the basis of simultaneous analysis of the three data samples. The branching ratio for annihilation 
into (pr] from the ^Si protonium state turns out to be about ten times smaller as compared to 
the one from the ^Pi state. An opposite trend was observed ||] for the (pn final state. 

We would like to thank the technical staff of the LEAR machine group for their support 
during the runs. We are grateful to J.Ellis, D.Kharzeev and V.Markushin for the stimulating 
discussions. 
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